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Short-range quark-quark correlations are introduced into the quark-meson coupling 
(QMC) model phenomenologically. We study the effect of the correlations on the structure 
' of the nucleon in dense nuclear matter. With the addition of correlations, the saturation 

curve for symmetric nuclear matter is much improved at high density. 

More than a decade ago, GuichonS' proposed a relativistic quark model for 
nuclear matter, where it consists of non-overlapping nucleon bags bound by the self- 
consistent exchange of scalar (a) and vector (o>) mesons in mean-field approximation 
(MFA). This model has been developed as the quark-meson couplint(QMC) model, 
and successfully applied to various phenomena in nuclear physics.™ Recently, the 
model has been further extended using relativistic constituent quark models. u 

So far, the use of the QMC model has been limited to the region of small to 
"^T* \ moderate densities because it has been assumed that the nucleon bags do not overlap. 

It is therefore of great interest to explore ways to extend the model to include 
short-range quark-quark interactions, which may occur when nucleon bags overlap 
at high density. In this report we will introduce these short-range correlations in a 
very simple way, and calculate their effect on the quark structure of the nucleon in 
medium. We refer to this model as the quark-meson coupling model with short-range 
correlations (QMCs). 

Let us consider uniformly distributed (iso-symmetric) nuclear matter with nu- 
clear density ps- At finite density the nucleon bags start to overlap with each other, 
and a quark in one nucleon may interact with quarks in other nucleons in the over- 
lapping region. Since the interaction between the quarks is short range, it may be 
reasonable to treat it in terms of contact interactions. An additional interaction 
term of the form, d n t ~ J2i^jJpq(^)^a''Pq(i)'ipq(j)r a ' t Pq(j), is then added to the orig- 
inal QMC lagrangian density.™ Here ip q (i) is a quark field in the i-th nucleon and 
r a stands for a combination of 7 matrices (with or without the isospin and color 
generators). In this report we shall consider only u and d quarks and simulate the 
short-range interaction using scalar- and vector-type couplings in MFA: r a = 1 and 
To- 

Next we consider the probability for the nucleon bags to overlap, using a simple 
geometrical approach. eP To measure the overlap of nucleons we treat them as the 
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MIT bag with the radius R. For any two nucleons the overlapping volume (measured 
in unitsnf the nucleon volume V/v( = 47i\R 3 /3)) is V ov (y) = 1 — 3y/4 + y 3 /16 (for 
y < 2),S and V ov vanishes beyond y = 2, where y = s/R with s a distance between 
the two centers of the nucleon bags. 

In a nucleus with A nucleons, any given nucleon may overlap with (^4 — 1) oth- 
ers. If the nucleons are distributed according to some two-particle density function 
P2(fi, r 2 ), then the overlapping volume per nucleon (in units of Vn) is 

V A = (A-l)f J dr 1 dr 2 p 2 (r 1 ,r 2 )V ov (\r 1 -r2\/R), (1) 

where the function p 2 may be given as the product of single-particle densities pi(r\) 
(normalized to unity), pi(r 2 ) and the two-nucleon correlation function F(\r\ — r 2 \). 
In the limit A — > oo, we then find the overlapping volume as Voo = (pbVn) x (1 — £) 
in the uniform matter, where 

Z = -LjdsH(s)V ov (s/R), (2) 

with H(s) = 1 — F(s). Here £ describes the effect of two-nucleon correlation in 
the overlapping part. If there were no correlation between two nucleons (namely, 
F(s) = 1 or £ = 0) Voo were 0.32 (when R = 0.8 fm and p = 0.15 fur 3 , the 
saturation density of nuclear matter) . 

For nucleons separated less than ~1 fm, some modification of the two-nucleon 
density is expected and this is described by the correlation function F(s). In princi- 
ple, one could calculate the correlation function within QMC self-consistently, as in 
the relativistic Brueckner-Bethe-Goldstone formalism etc.H* It is, however, not easy. 
In this exploratory study we shall use a phenomenological form of the correlation 
function. 

We use a convenient form of the correlation function proposed by Miller and 
Spencer 0) - so-called the Miller-Spencer correlation function. We then find 

H(s) = 1 - F(s) = l-(l- \g 2 p (y^j (1 + /( S )) 2 , (3) 

where g p describes the Pauli (exchange) correlations and is given by g p (y) = (3/y) ji (y), 
where y = sHf with the Fermi momentum kp. The function f(s) describes the strong 
repulsion in the core part and is parametrized as f(s) = — (1 — /3s 2 ) exp(— as 2 ), where 
a and j3 are parameters. EP 

The Bethe-Goldstone theory places a restriction on the overall size of f(s). The 
nuclear matter density times integral of the square of /(s) is known as the "wound" 
integral, k: 

k = Pb J dsf(s) 2 . (4) 

This quantity, which is the convergence parameter of the hole- line expansion, has a 
typical value of 0.12 around pb = Po-w A next condition on /(s) is obtained by the 
consistency requirement Q) 

|^(l-i 5 2 ( y ))(2/(s) + /(s) 2 ) = 0, (5) 
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to ensure the normalization on p2- Another quantity of interest is the correlation 
length l c defined bylZP 

/•oo 

£ c = - ds(2f(s) + f(s) 2 ). (6) 





We shall use the conditions (|5|) and @ to determine the parameters a and (5 and 
check the value of "wound" integral k numerically at the same time. 

Now it is natural to choose the probability of overlap to be proportional to V^,. 
(We here ignore higher-configurations like 3-body correlation etc.) In MFA the Dirac 
equation for a quark field in a nucleon bag is then given by 

[* 7 • d - K - g> - fsi^M) - CsSw + f2(4i) q }ho}i> q = o, (7) 

where m q is the bare quark mass, a and u are the mean-field values of the a and u 
mesons and g% and ifr are, respectively, the a- and u-quark coupling constants in the 
usual QMC model. cP The new coupling constants have been introduced for the 

scalar (vector)-type short-range correlations, and are given by f^,-, = (f%, ^/M 2 )VJx, 
(M = 939 MeV, the free nucleon mass). Note that since the coupling constants 
have dimension of (energy) -2 we introduce dimensionless coupling constants f^ v y 

In Eq.(|^), (ipqifiq) and (iplipq) are, respectively, the average values of the quark scalar 
density and quark density with respect to the nuclear ground state, which are ap- 
proximately^ given by the values at the center of the nucleon in local density approx- 
imation. [I'B 

Now we can solve the Dirac equation Eq.(|7]), as in the usual QMC, with the 
effective quark mass 



m 



i 



m q -(g$a + fl(^ q )), (8) 

instead of the bare quark mass. The Lorentz vector interaction shifts the nucleon 
energy in the medium |) 



e(fe) = y/M** + k* + 3(g«u> + f«(4%)), (9) 

where M* is the effective nucleon mass, which is given by the usual bag energy, ft 
The total energy per nucleon at density ps is then expressed as 



E, 



tot 



A rkF 1 

, , 3 / dWM** + fe 2 + 3(fl£w + f«(4<ip q )) + - — (mla 2 - m 2 ^ 2 ), 

(27T) d /9 B J H 2p B 

(10) 

where m a and m w are respectively the a and u meson masses. The u field created 
by the uniformly distributed nucleons is determined by baryon number conservation: 
u = Sg^pB/m^ = guPB/m^ (where g u = 3g%), while the a field is given by the 
thermodynamic condition: (dEtat/du) = 0. This gives the self-consistency condition 
(SCC) for the a fieldi 



a 



fdM*\ f k F M* , x 

dk Ttt^ , ( U ) 



(2vr) 3 m 2 V da ) J V 'M* 2 + fc 2 ' 



where (dM*/da) = -3g$S N (a) = -g a C N {a). Here g a = 3&S N (0) and C N {a) = 
Sn(&)/ Sn(Q), with the quark scalar charge defined by £V(cr) = J N dr V'gVv I n 
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Fig. 1. ^ as a function of ps- The solid (dot- 
ted) curve is for £ c = 0.75 (0.7) fm. 
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Fig. 2. Energy per nucleon for symmetric nu- 
clear matter. The dashed .curve gives the 
result of the original QMC. The solid and 
dot-dashed curves are for £ c = 0.75 and 0.7 
fm in QMCs, respectively. The region en- 
closed with the dotted curves is the empir- 
ical one. US' 



actual calculations, the quark density (iptipq) in the total energy may be replaced 
by 3/9b, and the quark scalar density, contributing to the effective quark mass, is 
approximately given as (ipqi/jq) = (rn^/g a )a because of the SCC.crEP 

Now we are in a position to show our main results. We choose m q = 5 MeV and 
the bag radius of the nucleon in free space Rq to be 0.8 fm. The bag constant B 
and usual parameter z, which accounts for the center of mass correction and gluon 
fluctuations, in the bag model are determined to fit the free nucleon mass with Rq 
= 0.8 fm - we find B 1 ^ = 170.0 MeV and z = 3.295.i The coupling constants g a 
and g u are determined so as to reproduce the average binding energy of symmetric 
nuclear matter (—15.7 MeV) at the saturation density po- 

In Fig. 1, we show £ in Eq.((2|). As an example, we set f§ = 40,/$ = 8 to fit 
the (observed) energy per nucleon at high pb (see Fig. 2). If we fix i c = 0.75 (0.7) 
fm, the conditions © and © give a = 1.05 (1.20) fm~ 2 , = 0.617 (0.708) fm~ 2 
and k = 0.12 (0.10) at p . We then find that gl = 50.09 (48.11), gl = 37.22 (33.80) 
and K = 360 (370) MeV. As shown in Fig. 1, £ gradually decreases and it becomes 
almost constant at high density. 

In Fig. 2, we present the saturation curve for symmetric nuclear matter. We 
can see that QMCs can provide the energy per nucleon, which lies in the empirical 
region (enclosed with the dotted curves). I12P It may imply the importance of the 
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Fig. 3. Scalar mean-field values in nuclear Fi §- 4 Ratio of the in-medium rms radius of 
medium. The curves are labeled as in the nucleon < (calculated using the quark 

p. 2 wave function) to that in free space (r q = 

0.5824 fm). The curves are labeled as in 

Fig. 2. 



short-range correlations at high density. If we choose a large value of f$ the energy 
becomes larger at high ps and it also makes the effective nucleon mass large. In 
QMCs we find M* /M ~ 0.84 at p$. We also show the scalar mean-field values at 
finite density in Fig. 3. In the original QMC the strength of the scalar field goes 
up to about 500 MeV at ps ~ 4poj while in QMCs it is much reduced at high ps 
because of the short-range (repulsive) interaction. 

Turning next to the size of the nucleon itself, as measured by the root-mean- 
square (rms) radius of the quark wave function, we see in Fig. 4 that the short-range 
correlations give a little enhancement. This effect is, however, not strong within the 
present parameter set, and the increase of the size around po lies in the upper-limit 
value analysed by the electron scattering experiment. HIP 

We summarize the role of scalar- and vector-type short-range correlations. The 
scalar-type correlation modifies the effective quark mass in a nuclear medium as (see 
Eq.®) 



m q 



fq m 2 



M- (12) 



In Eq,(|r2|) the mean- field part {g a <J term) provides a density dependence on the 
quark mass of order of O(pb) at low density, while the correlation part gives higher 
order contributions. This leads to the reduction of the scalar_mean-field value in 
matter, which may affect the in-medium nucleon properties. 13 '113) 

The vector-type correlation as well as the u meson shifts the total energy. The 
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oj and the correlation contribute to the energy as (see Eq.(|9|)) 



2 



3(^ + /^v ? )) = -f! 



m 2 , 



f q m 2 

1 + 9 kw- VM 



9B- (13) 



Since Voo depends on the density the correlation again gives the contribution of 
higher order in pg. It thus enhances the energy at high density. If we introduce 
the effect of higher-configurations like 3-body, 4-body etc., such correlations may 
have dependences of higher power of ps and play an important role not only at 
high density but also near the region of the phase transition to quark-gluon plasma. 
Those terms may correspond to higher order terms appearing in the chiral effective 
lagrangian for nuclear matter. [3 

In conclusion, we have studied the effect of short-range quark-quark correlations 
associated with nucleon overlap in MFA. We have found that a repulsive vector-type 
correlation makes the saturation curve close to the empirical one at high density. 
Furthermore, we have shown that the scalar- and vector- type correlations consider- 
ably modify the properties of nuclear matter at high density. While our inclusion of 
the correlations has been based on a quite simple, geometrical consideration, in the 
future we would hope to formulate the problem in a more sophisticated, dynamical 



way. 



Acknowledgements 



The authors greatfully thank AW. Thomas and G.A. Miller for valuable discus- 
sions. This work was supported by the Australian Research Council. 



[1] P. A.M. Guichon, Phys. Lett. B200 (1988), 235; See also, K. Saito and A.W. Thomas, 

Phys. Lett. B327 (1994), 9. 
[2] P. A.M. Guichon, K. Saito, E. Rodionov and A.W. Thomas, Nucl. Phys. A601 (1996), 349; 

K. Saito, K. Tsushima and A.W. Thomas, Nucl. Phys. A609 (1996), 339; Phys. Rev. C55 

(1997), 2637. 

[3] H. Shen and H. Toki, Phys. Rev. C61 f200 0L 045205: K. Sa ito. |tiucl-th/0010035[ 



[4] K. Saito, K. Ts ushima and A.W . Thomas, |nucl-th/9901084| (unpublished); I. Zakout and 
H.R. Jaqaman, |nucl-th/000303l| . 



[5] F.E. Close, R.L. Jaffe, R.G. Roberts and G.G. Ross, Phys. Rev. D31 (1985), 1004. 
[6] C.J. Horowitz and B.D. Serot, Nucl. Phys. A464 (1987), 613; G.A. Crawford and G.A. 

Miller, Phys. Rev. C36 (1987), 2569. 
[7] G.A. Miller and J.E. Spencer, Ann. of Phys. 100 (1976), 562. 
[8] B.D. Day, Rev. Mod. Phys. 39 (1967), 719. 
[9] P.G. Blunden and G.A. Miller, Phys. Rev. C54 (1996), 359. 

[10] M. Sano, M. Gyulassy, M. Wakai and Y. Kitazoe, Phys. Lett. B156 (1985), 27. 

[11] I. Sick, in: Proc. Int. Conf. on Weak and Electromagnetic Interactions in Nuclei, (Springer- 

Verlag, Berlin, 1986), p.415; D.H. Lu, A.W. Thomas, K. Tsushima, A.G. Williams and K. 

Saito, Phys. Lett. B417 (1998), 217. 
[12] K. Saito and A.W. Thomas, Phys. Rev. C51 (1995), 2757; Nucl. Phys. A574 (1994), 659. 
[13] G. Gelmini and B. Ritzi, Phys. Lett. B357 (1995), 431. 



